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We construct a new functional for the single particle Green's function, which is a variant of the 
^\ ' standard Baym Kadanoff functional. The stability of the stationary solutions to the new functional 

ON . is directly related to aspects of the irreducible particle hole interaction through the Bethe Salpeter 

ON 1 equation. A startling aspect of this functional is that it allows a simple and rigorous derivation of 

both the standard and extended dynamical mean field (DMFT) equations as stationary conditions. 
Though the DMFT equations were formerly obtained only in the limit of infinite lattice coordination, 
the new functional described in the work, presents a way of directly extending DMFT to finite 
dimensional systems, both on a lattice and in a continuum. Instabilities of the stationary solution 
at the bifurcation point of the functional, signal the appearance of a zero mode at the Mott transition 
which then couples to physical quantities resulting in divergences at the transition. 
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D ■ I. INTRODUCTION 

b 

Functional methods have long been used as a tool to study interacting systems. These involve techniques ranging 
i_J , from simple weak coupling perturbative methods to strong coupling expansions. Often, since it is not possible to 
sum the entire series of diagrams one takes recourse to various approximation schemes in the hope that they will 
capture some of the essential physics. In addition, functional methods, have a long history of systematizing these 
approximations. A direct way of testing the validity of any approximation is through the effective action approach to 
interacting systems. Typically, the effective action is obtained in terms of some relevant variable, for e.g. the density, 
magnetization, the order parameter, the single particle Greens function etc., as a perturbative series in the interaction. 
A well known and frequently used functional is the Baym Kadanoff functional r&fc (j^)El, which is the effective action 
l— "™ '. for the single particle Green's function G. Tf,k gives the free energy of the system at its stationary point which are 
I ' known to be saddle points and not extrema. In order to use variational approximations on these functionals, the 
functional should be convex. This is the case when one performs Legendre transformations with respect to time or 
£f) • frequency independent quantities as in density functional theoryn. Very recently, a functional of the Green's function 
! whose minimum gives the local Green's function of the Hubbard model on the Bethe lattice was investigated by one 
of us. This functional has several desirable properties, in particular it is a true extremum at the physical local Greens 
functionE. This motivates us to search for generalizations of the functional of Ref. [|to finite dimensions. 

In this paper, we present a new functional T new for the single particle Green's function, for which we can show that 
its saddle point is in fact an extremum, in the case that the irreducible vertex function in the particle hole channel has 
a definite sign (see below for a precise definition). This functional has many other useful properties primarily in the 
context of dynamical mean field methods which can now be naturally viewed as an approximation to this functional. 
The paper is organized as follows: first we introduce the Baym Kadanoff functional and our proposed functional 
and obtain the conditions under which the saddle points of this functional are extrema. We then proceed to show 
that the dynamical mean field theory ( DMFT) equations are a simple outcome of certain simple approximations 
on our functional, hence providing some justification for direct extensions of dynamical mean field theory to finite 
dimensions. We also discuss some potential applications pertaining to the Mott transition and the problem of adapting 
tight binding basis for dynamical mean field calculations. 
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II. THE BAYM KADANOFF FUNCTIONAL 



Consider the partition function Z (or equivalently the free energy W) of a system of electrons moving in a crystal 
potential V c (x) with Coulomb interactions V 

Z = exp[-/3W] = J D[tp^] exp[- J dxtP+(x)[d T - ^- + V c (x)]tp(x) 

- i / dxdx'ip + {x)i) + {x')V{x-x')ili{x')i){x)] (I) 
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Here x = (x, r) denotes the space- imaginary time coordinates. Our formulation is fairly general, it applies to both 
lattice models or continuum models. In the lattice case, the space coordinate x labels lattice sites and the integral 
symbol should be interpreted as a summation symbol over the lattice sites. We have omitted additional quantum 
numbers that can be carried by the field ip, such as spin and also orbital quantum numbers, in order not to clutter 
the notation. 

The effective action for the single Green's function G i.e, Baym Kadanoff functional, is defined as the Legendre 
transform of W with respect to the interacting single particle electron Green's function G(x,y) i.e., 

T bk [G] = W[J] + JG (2) 

where, W[J] denotes the free energy in the presence of an external source field J coupled to the electron Green's 
function. Using the condition J = STbk/SG to eliminate J in (||), we obtain the Baym Kadanoff functional 



bk 



[TrlogG - Tr[G 1 - G'^G + $[G\] (3) 



where, Go is the non-interacting Green's function and <&[G] is the sum of all two particle irreducible diagrams con- 
structed using the interaction and G. Stationarity yields the Dyson equation 

G- 1 = Go" 1 - E[G] (4) 

where T,[x,y) — 6$/6G(y,x)). At the saddle point (Q), Tbk is just the free energy of the system. An effective action 
approach is most useful when the saddle points are extrema. This is related to the stability of stationary solutions, 
which demands that leading order fluctuations around the stationary solution (0), increase the energy. To obtain a 
formal criterion for stability, we first consider the fluctuations around the stationary solution 

6T = T[G + 6G] - T[G] = f f SG(y, x) *T -6G(v, u) (5) 

J J 6G(x,y)6G{u,v) 

where T is any functional. In the functional treatment, it is essential to note that any quantity 0(x, y) is defined as 
0{x, y) = (a;|0|y) and hence should be viewed as a matrix in x and y. 

To examine the stability of a saddle point of the functional it is convenient to introduce a matrix G which acts on 
the Matsubara frequencies without affecting the spatial coordinates i.e., CG(ioj n , x, y) = G(—iw n ,x,y). Using the 
following property of Green's functions G(iw„,x, y)* = G(— iui n , x, y, which we assume to hold for the trial Green's 
functions as well, we can rewrite (||) as 

ST = y^6G(iu n )*x(iuj n , iu m )6G(iu) m ) (6) 

The matrix \ = G\ is defined by \ — S 2 T/SG 2 . While the frequency indices have been made explicit, matrix 
multiplication in the real or momentum space indices is implied. For the saddle point to be an extremum (maximum 
or minimum) the form defined by ffl) has to be negative (positive). For 



S 2 T bk 



5G{x,y)5G{i 



= -G-\y, u)G-\v, x) + TP h (xy- uv) (7) 



Here T ph (xy;uv) = (xy\T ph \vu) — 5 2 &/5G(x,y)8G(u,v), is the irreducible vertex which describes the irreducible 
particle hole interaction in the two particle channel and should be viewed as a tensor in the indices x,y,u,v. The 
analog of ([7]) in momentum space is 

Xbk(Pl,P2,P3,Pi) = -G~ 1 (p 2 ,P3)G~ 1 {p4,Pl) + T ph (p 1 ,P2,p 3 ,p 4 ) (8) 

where, p — (p, iu> p ) with ui p being the Matsubara frequencies. For T = T^fc, using (0) in (|^), we see that even in the non- 
interacting case, where $ is zero, the first term in (R) makes unbounded from below. To further illustrate this lack 
of stability in Tbk, we consider the simple example of an electron gas with Coulomb interactions V(x) = V(x)8(t). Let 
V(q) be the Fourier transform of U(x). In the Hartree approximation, $ is frequency independent and and rewriting 
the spin and momentum indices explicitly, 

F ph (pi,P2,P3,P4) = 5{pi +p 3 -p 2 -p4)U(p 3 -p 4 ) (9) 

Substituting this in (@) 
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Xbk(Pi,P2,Pa,Pi) = [-6(p2 -Pa)$(pi -Pi)G 1 (p 2 )G 1 (p 1 ) + S(p 1 + p 3 - p 2 - Pi)V(p 3 - p 4 )] (10) 



where, G~ 1 (p) = iui p — e p and e p is the dispersion of the particles shifted by the Hartree self energy term NV(0), 
with N being the total number of particles. 

We see from ( |l0| ) that Xbk, in the sector of zero total energy and momentum, Xbk decouples into 2x2 blocks 
in the non interacting limit, with zero diagonal matrix elements and equal off diagonal matrix elements ui n 2 + e 2 . 
Consequently, (ph has positive and negative eigenvalues even in the non interacting case. The Hartree term adds the 
constant NV(0) to all entries of Xbk ■ This results in the stationary solution remaining a saddle point rather than 
an extremum for small values of the interaction. Moreover, ( |Io|) shows that attractive interactions tend to further 
destabilize Tbk- 



III. A NEW FUNCTIONAL 



To construct a functional with a different stability matrix, we modify as follows: 

T new = T bk - Tr log(l - JG) - TrJG (11) 

In simple terms, we have changed the value of the functional away from the saddle point, by adding certain source 
dependent terms which vanish quadratically as the source is turned off. Note that though T new is closely related to 
Tbk, it is not a simple Legendre transform of the free energy in the presence of a source coupled to a simple composite 
operator. Using the relation J = 8T bk /8G = G -1 - G^ 1 + 6&/5G in ( pi] ) we get 

r nelu = -Trlog[G 1 -^]-TrG^ + 0> (12) 

It is straightforward to check that though T new and yield (^) and the same free energy at stationarity, their higher 
order derivatives evaluated at stationarity are different. This has profound implications for the stability. 

The stability of a stationary solution is again described by (j^, |^) , with T replaced by T new . In the non- interacting 
case, T new is trivially a constant. The second derivative of T new at stationarity is 

' H " ' = -T ph (xy]uv) + [ G{a,m)r ph (mn- 1 xy)G(n,b)T ph {ba:uv) (13) 



5G(x 1 y)5G(u,v) ' J abmn 

A diagrammatic representation of the above equation is shown in Figjl]. The integration over a variable x denotes 



u 

+ 



n 

FIG. 1. Diagrammatic representation of the (Il2 




In momentum space, ([13]) reads 



X ne W {hm 1 p,q) = -T ph {l,m,p 1 q)+ / G(fci, k 2 )T ph {k 2 , k 3 , 1, m)G(fc 3 , k uP , q) (14) 

J fei ,k 2 ,k 3 ,k4 

where Xnew is the Fourier transform of the second derivative of T new . In the translationally invariant case, defining 
the momentum transfer Q = I — m = q — p ( |l4] ) simplifies to (Fig. 2) 

Xne W {l,p; Q) = [-T ph il,p; Q) + J G(k - Q)r ph (l, k- Q)G{k + Q)r*> h (k,p; Q)\ (15) 

For the stability analysis we limit ourselves to consider the case of zero transfer frequency ujq = since it is natural 
to expect that the leading instability sets in first at zero frequency. This is known tOj-be the case of the two dimensional 
Hubbard model within the full parquet approximation as investigated by Bickersa. To derive the conditions under 
which the stationary point is an extremum, we analyze the form of Xnew , and then derive general stability conditions 

for it. The matrix T ph obeys the property r p,lt = CT ph C and so do (GG) and X new, implying that X n ew and CT ph 
are hermitian matrices. The particle hole kernel which enters in the Bethe Salpeter equation can be rewritten as 
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(1 — GGT ph ) — '}2 i ^i\<t>i){4 l i\ with all the eigenvalues positive as required by the stability requirements of the 
particle hole kernel. Moreover, x new = —'}2 i \iCT ph \(j) i ){(j) i \. If CT ph is positive definite (negative definite), (||) is 
negative (positive) definite and hence the stationary solution denotes a maximum (minimum). Applying the same 
argument to T^, we note that the functional is an extrcmum only if the Hermitian matrix CGG is positive or negative 
definite. Referring back to the previous section, we see that this is not the case even in the non-interacting limit where 
this operator has eigenvalues of both sign. To summarize, the instability of the stationary solutions of T new is linked 
to physical instabilities in the particle hole sector. Since the eigenvalues of the Bethe Salpeter equation are always 
positive when the solution is stable, instabilities are signaled by the vanishing of an eigenvalue for the first time. 
We now illustrate the stability analysis of T new in the Hartree approximation. Using (H) in (14), we see that 



X/r. 



,(l,m,p,q) =6(l-m+p-q)[ I G(k)G(k + 1 - m)V(l - m)V(m - 1) - V(l - m)] 

Jk 



(16) 



As before we will consider the case of zero frequency and momentum transfer. Summing over the internal frequencies 
iu>k in (|l^) we obtain 



,(l,p;Q) = F(Q)[-l + F(Q) 



Kk-Q)-n(k+Q)) 
k ( e (k-§)- e (k+§)) J 



(17) 



where, n(k) is the Fermi occupation function with the Hartree shifted chemical potential. 
Taking the Q — > limit, 



Xn 



,(1,p;0) = -K(0)[1-V(0) JJ^] 



Note that Xnew viewed as a matrix in Matsubara frequencies has constant entries, and so does Xr. 
of this matrix are zero with the exception of one single eigenvalue proportional to 

-V(0)[l + V(fl)N(ji)] 



(18) 

■j . All the eigenvalues 
(19) 



where N(p) is the density of states corresponding to the Hartree shifted chemical potential. From (19), we see that 
the Hartree approximation is again unstable for repulsive interactions. However, it is stable for attractive interactions 
V — > — V, and becomes unstable when V(0)N(ff) — > 1 triggering an instability in the particle-hole sector. This is 
the well known phase separation instability of fermions with attractive interactions. To investigate superconducting 
instabilities one would have to introduce anomalous Green's functions, and extend the recent work of Janiscl who 
considered the Baym Kadanoff functional, to our functional. A more comprehensive discussion of the stability of T new 
and the relative stability of T new and is left for future work. 
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FIG. 2. Diagrammatic representation of (Jl 5[) in momentum space 



IV. DYNAMICAL MEAN FIELD THEORY AS AN APPROXIMATION 



We now focus on the various applications of the new functional introduced in this paper. Despite the simple 
structure of both T{,k and T new , it is often impractical and intractable to sum over the infinite series of diagrams that 
contribute to 3>[G]. Rather, various approximation schemes are used to obtain (T new ) and hence G. The validity 
of these approximations depend on whether (7|) (or ( |l3| ) for T new ) are positive or not. Typical approximations range 
from Hartree Fock and RPA methods at weak coupling to atomic limits for strong coupling. Recent years have seen 
the development of Dynamical mean field theory (DMFT) which is one of the very few methods available to study 
correlated electron systems in the weak, strong and intermediate coupling regimes. Though this method is exact only 
in the limit of infinite dimensions, it has been very successful in describing the physics of realistic three dimensional 
systems, like the transition metal oxides. DMFT is formulated on a lattice and is characterized by self- consistency 
conditions (henceforth referred to as the DMFT equations) on the full local Green's functions. The form of these 
depend on the non-interacting density of states of the lattice considered. These equations were derived using the cavity 
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method or local perturbation theory in infinite dimensions!! Consequently, issues regarding the fluctuations around 
the DMFT, extension to finite dimensional lattices and systems in the continuum, have been difficult to analyze. 

Some of the above issues can be resolved if DMFT is viewed as an approximation to T new . Here, for the first time, 
we show that the DMFT equations can be derived directly from T new , by making the ansatz that G is completely 
local i.e., G(r, r', iu) n ) = G(r, r, iuj n )S(r — r'). Notice that our approach is fairly general. If we have in mind a lattice 
model the delta function <5(r — r') should be understood as a Kronecker delta function, while if we have a multiband 
model in mind Go and G should be viewed as matrices in the band and spin indices. Using the above ansatz, the 
stationarity condition for T new yields 

G(r,r,zc„) = [G - 1 -^|]- 1 (r,r) (20) 

where the second term within the square brackets is a diagonal matrix with entries 8<f>/8G = <5<!>/<5G(r, r, iui n ). The 
above equation should be viewed as a matrix in the spatial coordinates. Interestingly, ( p0| ) are exactly the DMFT 
equations postulated in Ref. |?] for non-translationally invariant systems. 
The local ansatz applied to yields the following equation 

G-\r,r,iu n ) = GoV.r.iu;,,) - . (21) 

ou(r, r) 



The crucial point is that (|2lj) and ( J2C| ) become inequivalent in the local G approximation. This becomes transparent in 
momentum space. Consider a translationally invariant system with Gg -1 = iuj n — tk where e& is the energy dispersion 
of the non-interacting system. Rewriting ( p0| ) in momentum space, 

- / T, Pie)de s, ; (22) 
where, G(r, r, i m n ) = G(iu> n ) and p(e) is the non-interacting density of states. On the other hand, using the above 



form of G 1 in (21) clearly leads to unphysical conditions. Replacing p in (|22| ) by a lattice density of states in the limit 



of large dimensions yields the standard DMFT equations. The potential applications need to be explored further. 



V. BASIS OPTIMIZATION CRITERIA IN DMFT CALCULATIONS 

Model Hamiltonians presuppose a choice of tight binding basis and hence model parameters. In a previous paperi, 
we argued that starting from the full Hamiltonian we can construct a functional of the local Green's function which 
yields the true local Green's function of the problem at stationarity. There, DMFT viewed as a theory for the local 
Green's function, is an exact theory but the functional which needs to be constructed is not known explicitly (beyond 
a general diagrammatic series discussed in refl) . Here we adopt a different approach and would like to view dynamical 
mean field theory as an approximation to a simple functional. When the saddle point fluctuations @ are positive, 
we can gauge the quality of different approximations. 

We first notice as in Ref. || that different choices of tight binding orbitals <p a where a denotes the band index, can 
be used to represent the exact one particle Green's function 

G(r,r',iuj n )= 0a(r-Ri)G Q/3 (Ri,R 2 ,^„)<^(r'-R 2 ) (23) 

Fti ,R.2 

We can then insert the following approximate trial Green's function into our functional 

G(r,r',^„)= ^ Q (r-R 1 )G Q/3 (Ri,R 2 ,^Mr'-R 2 )<5(Ri,R 2 ) (24) 

R.1 ,R,2 

compute the trial value of T new and optimize with respect to the G a p first by solving the DMFT equations which 
expresses the local lattice self energy 

£ Q/3 (Ri, Ri, iu n ) = 77; — 75 p — ; — r (25) 
dG Q . 1 a(Ri,Ri,zw n ) 

in terms of G Q| g (Ri, Ri, ioJ n ). Another choice of orbitals would lead to a different trial function (p4|), so a functional 
perspective offers the possibility of optimizing the basis of orbitals specifically for DMFT calculations. 
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VI. DERIVATION OF EXTENDED DYNAMICAL MEAN FIELD EQUATIONS 



Ordinary dynamical mean field theory treats longer range interactions at the level of the Hartree approximation. 
An extended dynamical mean field approach has been suggested! to partially remove this shortcoming. Recent work 
on the effects of long range Coulomb interactions within DMFT have shown that they can substantially modify 
the behavior near a Mott transition!^, changing the nature of the transition from a continuous one for short range 
interactions, to a discontinuous one. In the light of these works and given the relevance of long range interactions, we 



present a functional derivation of extended DMFT equations along the lines of Sec. IV. 

We consider the system (Q) for electrons with some non-local interaction. Decoupling the repulsive quartic interac- 
tion term in (Q) using a Hubbard Stratanovich field <f>, we get 

Z = J D[tp^(j)] exp[- J dx^ + (x)[d T - + V x (x)]ip(x) 

-~ / dxdy^x)^ 1 (x - y)<p{y) ~ i J dx<f>(x)^ (x)ip(x) (26) 

If we chose Hq(x, y) — V{x — y) to be instantaneous in time, then upon integrating out the the <f> field we recover the 
interacting electron system of (|l|). This decoupling can be extended to attractive interactions between the electrons 
by retaining the same repulsive interaction between the phonons (so as to obtain a positive dispersion) and replacing 
the factor i in (|26|) by unity. 

Physically, 4> can be interpreted as a phonon field with a particular dispersion determined by IIo. In fact, we can start 
with (p6| ) as a general action describing interacting electrons and phonons, the limit of static Coulomb interactions 
being obtained as a limiting case. In momentum space, the term quadratic in the phononic field is (^(q)!!^ 1 (q)<p(— q) 
with IIo being the Fourier transform of the interaction. We now have a system of interacting fermionic and bosonic 
fields. Similar to Secjll] , we calculate the free energy W[J,K] of the electron-phonon system, where J and K are the 
source fields for the electron Green's function G and the phonon Green's function II respectively. Note that in the 
present case, the phonon Green's function is nothing but the full density density correlator or two particle Green's 
function of the electrons. The effective action functional i.e., the equivalent of the Baym Kadanoff functional for such 
a coupled system is now given by the double Legendre transform of W[J, K] 

T[G,n] =TrlogG-Tr[G ( 7 1 -G^^G- ^Tr log n+iTrn - 1 n + $[G,n] (27) 

<fr is the sum of all two particle irreducible diagrams constructed using the electron-phonon vertex. The source fields 
are fixed by the following conditions: 

J = G- 1 - Go" 1 + g 

k = rr 1 - n 1 - 2^ (28) 



The saddle point conditions on the above functional are given by (|28|) with J and K set equal to zero. The stability 
condition for this functional is given by 

/ f^X^x) — -6X 1 (u,v)>0 (29) 

J J Tl oX i {x,y)8X J {y,u) 

where Xi refer to G and II. Following (|ll|) we introduce a slightly modified functional 

T new [G,U] = r - Trlog(l - JG) - TrJG + -[Tr log(l - ATI) + TrKU] (30) 



Substituting (ES) in (BOh, we get 



,[G,H] = -Tr log ( G^ 1 - — J - TrG— + -Trlog f IL7 +2- - TrTI— + $[G,H] (31) 



These extra terms tend to stabilize the functional by adding positive terms to T. The saddle point equations of ( |3l| ) 
are 



G 



SG 

sn 



"0 8G 



-G 



G n 



8<S> 
SG 



G 



5 2 4> 
SG5G 

<5 2 $ 

susg 



IE 



o (5<I> 

z «5n 



n- 



-n 



n 



6 2 <P 
SGSU 

5 2 $ 

susn 



(32) 
(33) 



The saddle point solutions are obtained by equating the above equations to zero. The relevant solution is G 1 = 
5&/5G and IT 1 = 1 + 25<$>/5D since © demands that when J = K = 0, <ST new /0: = (T/tSX. The second 
order derivatives are given by 



Go 1 



5 2 T r . 



5GSG 
S 2 T 

u L new 

smn 

S 2 T 

u L new 

SG5T1 



S 2 ® 



SG5G 5GSG 
5 2 <S> 8 2 $> 

2n— — n 



SG5G 
6 2 ® 



S 2 ® 8 2 ® 

2n— — n- 



susu susu 



6 2 3> 8 2 <S> 
n n 9TT 

SGSn SGSG 



SG5U SGSU 
q ^ q 5 2 ® 



5USU 

5 2 <f> 

snsn n 5Gsii ~ SGsn 



SG6U SGSG 
5 2 <S> <5 2 $ 



To obtain the extended mean field equations, we make the ansatz: G(r, r',iw„) 
n(r, r', iu n ) = n(r, r, iu> n )S(r — r'). Minimizing with respect to G and II we get 



(34) 
(35) 
(36) 

G(r, r, iuj n )S(r — r') and 



G(r,r,iu n ) = [G 1 

n(r,r,^„) = [n - x + 2 



<5$ 
~5G 



<5IT 



l (r,r) 
>,r) 



(37) 



As before, 5Q/5G and 6&/5H are local self energies and should be viewed as diagonal matrices in the space coordinates. 
In the translationally invariant case, the above reduce to the following equations 



G(iu! r , 



r/q 



U{iuj n ) = / — - 



f/q 



(q) + 25<i>/Sn 



(38) 
(39) 



These are precisely the extended dynamical mean field equations obtained in Ref. We therefore, see that all the 
dynamical mean field theories can be seen as consequences of local approximations made on the new functionals. An 
analysis of stability similar to the one leading to the conditions on (0) is left for future work. 



A. Mott transition 

Recently, the Mott (metal-insulator) transition in the infinite dimensional Hubbard model! was revisited from a 
Landau Ginzburg perspectives!. The central idea consisted of the introduction of a Landau functional whose saddle 
point gives the (local) Green's function of the problem. The Mott transition is then viewed as bifurcations of the 
stationary points of this functional as the parameters of the theory are varied. This approach was used to clarify various 
issues related to the zero temperature Mott transition in infinite dimensions, and to study the finite temperature Mott 
transition at half filling in the infinite dimensional Hubbard modeLJ. 

Here we suggest that the functional T new discussed in this paper, be viewed as a generalization of the Landau 
functional discussed in Ref. || to finite dimensions, and use it to extend some of the lessons learnt in infinite dimensions 
to finite dimensions. In particular, we elaborate on a connection between the existence of a Mott transition and the 
divergence of the. charge compressibility, conjectured on the basis of pioneering numerical work by Furukawa and 
Imadall3. TheyL3 discovered that the approach to the Mott transition in two dimensions as a function of doping 
is characterized by a divergent compressibility. By varying the hopping integrals, they demonstrated that this is a 
general phenomena characteristic of the approach to the Mott transition. Here we use the following argumentEj. The 
Mott transition is but a bifurcation of the stationary solutions of T new i.e. solutions of ST new [G]/SG = 0, as a control 
parameter a (which could stand for chemical potential , temperature or U ) is varied. For a bifurcation to take place 
the matrix \ new — &Gl ^ x y ^YG(u v) nas ^° ^ e non invertible since otherwise the equation: 



7 



4 °-(w)r(w). & ' w 

could be uniquely solved for 5G. In (fh)|), the subscript c denotes the quantities evaluated at criticality. We have seen in 
Sec.[il]| that when Xnew ceases to be invertible, an instability sets in the particle hole channel at zero momentum. This 
implies a divergent compressibility unless there exist matrix elements which cancel the contribution of the vanishing 
eigenvalue of the Bethe Salpeter kernel, which signals the instability. This argument directly relates the existence of 
a Mott transition to the divergence of the compressibility supporting the observation of Ref. [l^. We stress however 
that while our argument is fairly general and very appealing in the light of the numerical results of Ref. |l2|, it is based 
on assumptions which need to be verified on a case by case basis. In fact, even in infinite dimensions, the Landau 
functional at zero temperature has some very singular bifurcations in its phase diagram at which the compressibility 
does not diverge. This issue as well as the computation of matrix elements in specific bifurcations requires a problem 
specific analysis which at this stage can be performed rigorously only in the limit of infinite dimensionstS. 

Finally it would be useful to extend the analysis to spin dependent Green's functions. This would allow us to 
understand the effects of magnetism in many cases, examples being the doping driven Mott transition in two di- 
mensions, where one approaches a magnetically ordered state and the problem of of the antiferromagnetic metal to 
antiferromagnetic insulator transition in infinite dimensions that we have studied recentlyo. 



VII. CONCLUSION 



To summarize, we have constructed a new functional for the single particle Green's function and investigated its 
use as a variational free energy. We found that stability of the saddle point solutions was linked to the Bethe-Salpeter 
equation in the particle hole sector. From the perspective of applications, the use of V new allows us to derive DMFT 
and extended DMFT, directly in finite dimensions, for non translationally invariant geometries. Functional methods 
provide simple connections between one particle Green's functions and particle hole instabilities, and we use these 
ideas to discuss the compressibility near the Mott transition. Finally, we have shown that under some circumstances 
it is possible to understand whether the stationary points of T new are extrema or saddle points. We feel that this is a 
significant advantage, over the Baym-Kadanoff functional Tf,k, when used to generate non perturbative approximations 
using trial Green functions. There is a great deal of freedom, in adding sources to so as to convert its stationary 
points into extrema, and constructions along the lines of our paper, but with different additions of sources, so as to 
obtain more general stability conditions than the ones we were able to obtain are well worth pursuing. 
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